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We investigate the spectral statistics of a network model for a three dimensional layered quan- 
tum Hall system numerically. The scaling of the quantity Jo = |(s 2 ) is used to determine the 
critical exponent v for several interlayer coupling strengths. Furthermore, we determine the level 
spacing distribution P(s) as well as the spectral compressibility \ at criticality. We show that 
the tail of P(s) decays as exp(— ks) with k, = 1/(2%) and also numerically verify the equation 
X — (d — D2)/{2d), where D2 is the correlation dimension and d — 3 the spatial dimension. 

KEYWORDS: network model, layered system, level spacing distribution, spectral compressibility, scaling, quantun 
Hall systems, multifractality 
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Although the quantized Hall effect is usually attributed 
to two-dimensional (2D) systemsEP the possibility of an 
occurrence in three-dimensional (3DA,systems has been 
theoretically investigated early onBtP In addition- 
some precursory experimental findings in 3D systems! 
there have been two classes of quasi-3D systems that 
show the integer quantum Hall effect: molecular crys- 
tals? and multi-layer quantum wpJls formed by graded 
GaAS/AlGaAs heterostructures.LrB' Recently, the latter 
have attracted significant theoretical interest. The ex- 
istence of a metallic phase for a finite range of energies 
and a new c/wmL2D_metallic phase on the surface have 
been proposedBt9'ELlP In ref. ^ these 3D systems of 2D 
quantum Hall layers coupled by weak interlayer tunnel- 
ing have been investigated using a network model. The 
existence of three phases: insulator, metal and quantized 
Hall conductor accompanied by extended surface states 
was shown by performing numerical transfer matrix cal- 
culations within the network model. 

We will use the same kind of network model to in- 
vestigate spectral properties of the 3D layered system. 
In doing so we will apply the methods-used for the 2D 
Chalker-Coodington network modefiaEif which use the 
spectrum of the unitary network operator U to determine 
the statistics-joLthe eigenvalue spectrum of the underly- 
ing system Jl3-t3 This method has the advantage that 
after determining the critical energy E c one can use the 
entire spectrum of eigenphases of U(E C ) to investigate 
the critical spectral statistics. 

The paper is structured in the following way: After an 
introduction to the network model we will use finite size 
scaling of the function J (L,E) — ^(s 2 ), where s is the 
level spacing in units of the average level spacing A, to 
determine the critical point. At the critical point we are 
going to investigate the level spacing distribution func- 
tion P(s) and the number variance T,2(N) = ((n— (n)) 2 ) 
of an energy interval containing on average N — (n) lev- 
els. We will see that, as in the 2D case, the following 



equation proposed in r ef. holds: 



X 



2d' 



(1) 



where \ = limjv-Kx) lirai,-^ dEa(iV) / 'dN is the spectral 
compressibility, d the spatial dimension of the system 
and 77 the diffusion exponent, which isolated to the 
fractal correlation dimension D2 = d — r]S3 We are also 
going to show that the tail region of the level spacing 
distribution shows an exponential decay at criticality, i.e. 
P(s) oc exp(— ks) and that as proposed in ref. [19f 
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The network model for the 3D layered system devel- 
oped by Chalker and Dohmen (CD) is the natural gen- 
eralization of the 2D Chalker-Coddington (CC) network 
modelBE 2 ! It is a model for electrons in 3D moving in a 
smooth disorder potential in 2D layers with a strong per- 
pendicular magnetic field. In the classical limit, reached 
when the correlation length ly of the disorder potential 
V is much larger than the magnetic length i c , the motion 
of the electrons can be decoupled into three independent 
components. The first with the smallest time scale is the 
cyclotron motion about a guiding center. Secondly, we 
have oscillations parallel to the magnetic field on an in- 
termediate time scale and for long times we have a drift 
motion along the equipotential contours of V in the plane 
perpendicular to the magnetic field. 

The network model simulates the system in the follow- 
ing way (see also Figure I): The unidirectional channels 
in the 2D planes represent the equipotential contours of 
the potential in the plane. The random length of the 
contour is simulated by a random phase added to the 
electron at each channel link. The electrons travel along 
those lines until they reach a node of the network which 
represents a saddle point of the potential. At the sad- 
dle point they are scattered. The scattering amplitude 
T = ( I + exp ((E — u) I 'E t ) is determined by the electrons 
energy E and the saddle point energy u, where E t is the 
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Fig. 1. The 3D Chalkcr-Dohmcn network. It consists of layers 
of the 2D Chalker-Coddington model coupled by interlayer scat- 
tering matrices. The dashed double arrows mark the interlayer 
coupling nodes. 
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characteristic tunneling energy of the potential 
the following all the saddle point energie^wilLbe set to 
zero to avoid classical percolation effectsJiSEJcir In ad- 
dition to the scattering at the nodes at each link a scat- 
tering in the direction perpendicular to the plane takes 
place. This is accomplished by a scattering matrix of the 
form 



cos <p 
— sin < 



sin (, 

COS( 



(3) 



giving us a transmission rate of t = sin 0, where 4> is the 
same for all interlayer scatterings. 

A state or wave function of the network is defined by 
a normalized vector = (ipi, . . . , ip n ) whose elements 
ipi denote the complex amplitudes on the n — 2L 3 net- 
work channels, where L is the system size. Its time evo- 
lution is determined by the scattering matrices at the 
saddle points and the interlayer couplings which can be 
expressed by one unitary matrix U the so called network 
operator: 



*(t + 1) = J7*(r). 



(4) 



Actually, U can be written as the product of a matrix 
U2d built from the independent network operators of the 
respective layer systems and an interlayer scattering ma- 
trix S which takes care of the coupling of those systems, 
i.e. U = UzdS- Both matrices are sparse with only 2 
non-zero elements per column, which means that U is 
also sparse with 4 non-zero elements per column. U is 
energy dependent through the energy dependence of the 
scattering matrices in the layers. 

An energy eigenstate ^ e of the underlying electronic 
system is a solution of the stationarity equation: 



U{E)^! 



(5) 



which corresponds to an eigenstate of U(E) with eigen- 
value unity. Such an eigenstate will occur only at discrete 
values of E = E c which form the energy spectrum of the 
modeled system. However, it has been shown that the 
eigenphases u> a (E) of the unitary network operator U{E) 
determined by the equation 



U(E)* a 



(6) 



show the same statistics as eigenenergies E n close to 
Er3 Since the arguments in r ef.|l4|, given for the Chalker- 



Coddington model, hold for any network model that in- 
cludes an energy parameter we can use the so called 
quasienergies Lo a {E), which can be considered an exci- 
tation spectrum at energy E, to determine the spectral 
properties of the 3D system. 

Since the 3D network model exhibits three different 
phases, i.e. two phase transitions, our first task is to 
determine the critical points of those transitions. It 
has already been established in ref. |9| that at energy 
E Cl (t)(< 0) the network exhibits an insulator-metal (IM) 
transition and symmetrical to the E = axis at en- 
ergy E C2 (t) = —E Cl (t) a metal-quantum Hall conductor 
(MQH) transition for every t. The states in the quan- 
tum Hall conductor are all localized as in the insulator 
phase except for the fact that for hard wall boundary 
conditions there exist extended surface states separated 
from the localized bulk states. In our investigation we 
are using mostly periodic boundary conditions, but we 
will comment on the hard walled systems if results differ 
significantly. 

The width W = E C2 — E C1 of the metallic phase was 
shown to depend on the coupling t in the following way: 
W(t) oc i 1 /" 2 " 1 , where V2d ~ 2.3 is the critical exponent 
of the localization length in 2D quantum Hall Systems HP 
This knowledge was used to narrow the range for our 
search of the critical energies. 
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Fig. 2. The scaling variable Jo as a function of E for different 
system sizes and t = 0.2. At E rs 2.2 all the curves meet in a 
single point. The inset shows the one-parameter dependence of 
Jo on L/£(E). 



The critical energies are determined by finite size scal- 
ing. As a scaling variable the quantity Jo = \ (s 2 ) is used. 
This function reflects the transition of the spectral corre- 
lations from Gaussian unitary ensemble (CUE) to Pois- 
son statistics when we cross from the metallic phase into 
a localized regime. The values Jo = 0.590 in the metallic 
and Jo = 1 in the localized regime are well known from 
random matrix theory (RMT).o 

In Figs. 2 and 3 we can see the function Jo(E) at t = 
0.2 for different system sizes L and energies E > and 
E < 0, respectively. The curves for different L in both 
figures intersect in single points at E Cl ~ —2.2 and E C2 
2.2 with Jq p» 0.68 in both cases. These points mark 
the IM and MQH transitions. Assuming the validity of 
one-parameter scaling, i.e. Jq(E 7 L) — f(L/£(E)), and 
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Fig. 3. The scaling variable Jo as a function of E for different 
system sizes and t = 0.2. At E —2.2 all the curves meet in a 
single point. The inset shows the one-parameter dependence of 
Jo on L/$(E), 



the existence of the localization length £(E) that shows 
the power law behavior £(E) oc \E — -E C | _I/ , one can re- 
scale L so that the data of Jo as a function of L/£ for 
consecutive energies E overlap with each other. Within 
statistical accuracy all points fall onto one common curve 
independent of L and E. This can be seen in the insets of 
the two figures where Jq(E, L) is plotted against L/£(E). 
In order to calculate the critical exponent v we expand 
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Table I. The results for the 4-parameter fit of eq. (I 



pare it to the Poisson distribution (Pp( a ) = exp(— s)) 
in the localized and the GUE distribution (Pgue(s) = 
|t s 2 exp ( — § s 2 ) ) in the metallic regime known from RMT 
The critical P(s) shows characteristics of both of these 
distributions. In Fig. 4 we can see the critical P(s) for 
t = 0.2 and several system sizes. As expected they are 
independent of the system size and all fall on the same 
curve. They also obviously differ strongly from the Pois- 
son and GUE distributions. Those distributions were ob- 
tained from networks with periodic boundary conditions. 
The inset of Fig. 4 shows a comparison of the periodic 
boundary case with a critical distribution of a system 
with hard wall boundaries (HWB) in one direction. We 
can see that in the HWB case P(s) is shifted towards the 
Poisson distribution which means a different 0.7. 
This has also been observed in other 3D systemsJ 2 3 ) We 
also found a slightly increased value for E c which if con- 
firmed has to be investigated in later publications (see 
also refs. ||and|3C]). 
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Fig. 4. The critical level spacing distribution P(s) for t = 0.2 
and several system sizes compared to the Poisson and GUE dis- 
tributions. The inset shows the critical P(s) for two different 
boundary conditions. 



the finite-size scaling function around the critical point: 

J {E,L)^4 + B(E-E c )L 1 /». (7) 

We can now either determine v, Jq, E c and B by a four- 
parameter fit or we can take the energy derivative of 
eq. (M) , which gives us dJ (E, L)/dE oc L x l v and perform 
a two parameter linear fit to a double-log plot of the 
slope of J(E) as a function of L at the critical point 
E c . The results for the 4-parameter fit are shown in 
Table I for different values of the coupling t. The linear 
fit of the slope data yields essentially the same results 
as the 4-parameter fit (see Table II). The results for 
the exponent v are in agreernpatwith earlier numerical 
results for layered 3D systemsE-rtf [y — I.3I1JZ.Q— 15 and 
v = 1.45 ± 0.25) and 3D Anderson modelsfel'BEZP {u 
1 .4.'! ■ O.ilfi and v = 1.34±0.I0) with broken time reversal 
symmetry. The values for the critical energy also agree 
with the results of Chalker and Dohmen.EP 

Now that we know the critical energies E c (t) we can 
concentrate on the investigation of the spectral proper- 
ties at criticality. We start by looking at the critical 
level spacing distribution P(s). We are going to com- 



A closer investigation of the shape of the critical dis- 
tribution reveals that for small values of s it resembles 
the GUE case. The inset in Fig. 5 shows the small s 
behavior and we can see that P(s) oc s 2 for all couplings 
t. This reflects the level repulsion due to the multifrac- 
tal eigenstates which, like the homogeneously smeared 
out metallic eigenstates, cover the entire system. On the 
other hand, a look at the behavior for large values of s 
shows an exponential decay of the tail of P(s) (see Fig. 5) 
which is typical for the Poisson distribution and in con- 
trast to the Gaussian tail of the GUE case. Compared to 
the Poisson case, the critical distribution decays faster 
with P(s) oc exp(-Ks) where k = 1.76 ±0.04 for t = 0.2, 
which is the coupling where we got the most statistics. 
According to ref. |l^, which predicted this exponential 
decay, the factor n is connected to the spectral com- 
pressibility x Via eq. (g). This leads us to the next step 
in our investigation, the determination of x- 

In the right inset of fig. 5 we can see the level number 
variance E 2 whose slope gives us x- The results for the 
different coupling strengths all agree with each other and 
we again took the value for x from our data for t = 0.2. 
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Fig. 5. The tail of the critical P(s) for t = 0.2, 0.4, 0.5 compared 
to the GUE and Poisson case. The left inset shows the small s 
behavior of P(s) compared to the function s 2 . The right inset is 
the level number variance at criticality. 



We obtained: x = 0.286 ± 0.003 which gives us accord- 
ing to eq. (§)« = 1.75 ±0.01 which is in good agreement 
with the directly obtained value of k. These values for \ 
and k are also quite close to the ones obtained in ref. |3l] 
for a 3D Anderson model with time reversal symmetry. 
Furthermore, we can connect the spectral compressibility 
directly to the multifractal critical wave functions. This 
was done in ref. |l^ resulting in eq. ([!]). If we use x we 
obtain a value of ry = 1.72 ± 0.02 or D 2 = 1.28 ± 0.02. If 
we compare D 2 with direct results for the network model 
(1.3±0.ll23') or 3D Andeppn modfils without time rever- 
sal symmetry (1.7 ± 0.2j2j'w 1.5E-j) we see that eq. (|l|) 
holds within the errors. Although the result for D 2 taken 
from cq. |l], which fits very well to the direct result from 
ref. |3^, is a little smaller than those of other investiga- 
tions. 

In conclusion, the numerical investigation of the spec- 
tral statistics of a network model for a three dimen- 
sional layered quantum Hall system for different coupling 
strengths between the layers showed that the phase tran- 
sitions from metallic to localized regimes all have the 
same characteristics. The critical exponent v w 1.3 is 
in agreement with other results for similar systems. Fur- 
thermore, the relations (|l|) and (|J) were verified to a high 
degree of accuracy within the network model. This is an 
improvement to other 3D models (all of them with time 
reversal symmetry) where the two equations could also 
be verified but with much larger errorsjl3- ) This shows 
that the usage of network models, apart from their ad- 
vantages for the investigation of time evolution, gives us 
a more efficient way to investigate spectral properties. 

The author would like to thank Y. Ono for valuable 
discussions and the Deutsche Forschungsgemeinschaft (DFG) 
for their financial support. 
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Table II. The results for v obtained from the derivatives d Jo /dE 
at E c . 



